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Abstract 

Based on a first order gradient expansion a consistent transport equation is 
derived for a nonrelativistic system beyond the quasiparticle approximation, 
i.e. for a regime where the dynamically generated width of the states is allowed 
to be large. An exactly conserved quantity is identified which is interpreted 
as an effective particle number obtained by coarse graining. Using a test 
particle ansatz for this conserved quantity allows to rewrite the transport 
equation into equations of motion for test particles. The two-body collision 
terms are formulated in terms of the test particles which gain non-trivial 
renormalization factors due to the coarse graining process. 

PACS numbers: 24.10. Cn, 05.60.-k, 05.70.Ln 
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I. INTRODUCTION 



Since the pioneering works on non-equilibrium quantum field theory JI|-|J semiclassical 
transport theory has become a major tool to solve various problems of many-body physics. 
While former works have focused their attention more or less on the quasi-particle regime 
(see e.g. |5|-[TT| and references therein) and small corrections to it [ lT2|— p~3 1 the extension of the 
formalism to off-shell phenomena has become a topic of growing interest in the last few years. 
E.g. in heavy-ion collisions it has turned out that the collision rates usually are so high that 



an on-shell approximation seems to be inappropriate |T5Hl8|. In addition, the resonances 
excited during the reaction may have large decay width. Therefore, a representation of these 
states by stable particles may not be a proper approximation. 

In this article we shall set up a formalism which allows for a consistent transport the- 
oretical description of states with a continuous mass spectrum. For simplicity we restrict 
ourselves to the non-relativistic case. We will comment in the last section on possible gener- 
alizations to relativistic systems. We will derive a transport equation in first order gradient 
expansion which includes off-shell effects. We will discuss in detail how this transport equa- 
tion can be solved by a proper test particle ansatz. To account properly for the off-shell 
behavior the test particles are allowed to have an arbitrary energy not restricted by the 
mass-shell condition. While a test particle representation of a transport equation in the 
quasi-particle limit is straightforward we will see that this becomes a rather non-trivial issue 
if off-shell effects are included. Especially we have to worry about a quantity which is fully 
conserved by the transport equation. Only such a quantity can be represented by test par- 
ticles. Otherwise the test particle ansatz would yield a number of equations which is higher 
than the number of test particle coordinates, i.e. an overdetermined system of equations. 
We will isolate a quantity which is fully conserved by the transport equation and derive the 
corresponding test particle equations of motion. Finally we will comment on the form of the 
collision terms for these test particles. 

The article is organized as follows. In the next section we discuss in some detail for a 
simple example how a test particle ansatz works for a conserved quantity and fails for a 
non-conserved one. In Sec. |TT| equations of motion for the quantities of interest are derived 
from the underlying quantum field theoretical Kadanoff-Baym equations by performing a 
first order gradient expansion. This derivation is very similar to the one presented recently 
in |16| (second reference). Therefore, Sec. |T| might be seen as a review. It is included to 
keep the article self-contained and to put emphasis on various other aspects as compared 
to previous works. In Sec. [TV] new material is presented. Especially a quantity is isolated 
which is conserved by the transport equation derived in the preceeding section. For this 
quantity a test particle ansatz is presented in Sec. [V] and the equations of motion for the 
test particles are derived. Finally we evaluate the self-energies in the Born approximation 
using the test particle representation. We will summarize our results in Sec. [VT| , compare 
them to previous works and give an outlook on unresolved problems. 
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II. TEST PARTICLE ANSATZ 



Before turning to the proper transport equation for modes with a continuous mass spec- 
trum we would like to discuss how a test particle ansatz for the quantity of interest works 
and for which cases it makes sense. For this aim let us consider an equation for a quantity 
f(t,p) which is too complicated to solve it exactly. The only thing we have to know about 
this exact equation is that it conserves a particle number 

N(t) = Jdpf(t,p). (2.1) 

By an approximation scheme one obtains from the exact equation the following approximate 
one (which resembles a Vlasov equation): 

^ [(1 - K(t,p))f(t,p)} + F(t) ^f(t,p) = (2.2) 

with a force term F(t) and a renormalization k which both are not further specified. Ob- 
viously this approximate (transport) equation does not conserve the particle number N but 
instead the quantity 

N(t) = Jdp(l-K(t,p))f(t,p). (2.3) 

We now discuss two cases to solve ( |2.2|) , namely test particle ansatze for / and for / = 
(1 — re)/, respectively. In general, a test particle representation for any of the two quantities 
is given by 

/ 
or 

/ 

where L is a normalization constant, % numbers the M test particles and Pi(t) characterizes 
the trajectory of the test particle i. Obviously the momentum integral over Q2.4|) is conserved 
and given by M/L. Therefore, a test particle ansatz for / to solve ( |2.2j ) seems to be 
inappropriate since the momentum integral over / (which yields the particle number) is not 
conserved by ( |2.2| ). Nonetheless, it is instructive to figure out which equations one gets with 
a test particle ansatz for /. In this case one finds 



-i M 

7 £*(p-ft(*)) (2-4) 
L i=i 



= E |_^| i M, (p _ pi(i)) + 



^{p- Pi (t))\- (2-5) 



Note the appearance of the total derivative with respect to t in contrast to the partial 
derivative in (p.2|) . This is obtained by 



|[(1 - K(t,p)) S{p - ft)] = - 9 -^ L 5{p - ft) + (1 - K{t,p)) p{p - Pl 

(h;( '- p) Sip - ft) + (1 - K {t,p)) ^8ip - Pl ) (-ft) 



dt dp 
-j^- 1 - 8(p - Pi) - Pi dp 1 8(p-Pi) - (1 - K(t,Pi))pi—5(p-pi) 

(U;[ ' ' lh) 5ip - ft) - (1 - K (t, ft)) ft ^5ip - ft) (2.6) 



dt op 
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where we have used the identity 



, , . d , , . . d . dhiy) c/ . 

h{x) -q^°( x -y) = h[y) fa°( x ~y> — dy~ ^ x ~ y ' 



(2.7) 



To fulfill (|2.5|) one has to demand that the coefficient of the 5-function as well as the coeffi- 
cient of the momentum derivative of the 5-function both vanish, i.e. 



_ dK(t,Pi(t)) 
dt 

F(t)-[l-K(t, Pi (t))]pi(t) 








(2.8) 
(2.9) 



In this way we obtain two equations of motion for one test particle coordinate: 



Pi{t) 
Pi{t) 



dpi 
Fit) 



1 - K(t,Pi 



(2.10) 



(2.11) 



Obviously, this provides an overdetermined system which is simply due to the fact that a 
test particle ansatz has been made for a quantity which is not conserved by the equation 
one wants to solve. On the other hand, a test particle ansatz for / yields 



o = E 



F(t) 



K(t,Pi(t)) 



Pi(t) 



dp^^ 



Pi(t)) 



(2.12) 



Here the coefficient of the pure 5-function vanishes. This is a general property if a test 
particle ansatz for a conserved quantity is made. We only have to fulfill one equation: 



Pi(t) 



F(t) 



1 - nit, pi) 



(2.13) 



Here we observe a second interesting feature of a test particle ansatz. The last equation 
is identical to the second one of the equations obtained before. Both equations fl2.11|) and 
( |2.13| ) were obtained by demanding that the respective coefficient of the d p 5{p — pi) term has 
to vanish. Thus, the difference between a test particle ansatz for a conserved quantity and 
for a non-conserved quantity lies only in the appearance of a term with a pure 5-function 
(i.e. without derivatives acting on it). This provides an additional equation which causes 
the system of equations to be overdetermined for the case of a non-conserved quantity. The 
other equations do not differ for both cases. 

So far we were only concerned with a proper solution of the transport equation (|2.2|) and 
indeed we have presented one consistent way to solve it, namely by a test particle ansatz 
for /. Let us now recall that ( |2.2|) is only an approximate equation and that the exact 
equation in contrast to the transport equation exactly conserves /. Thus, one has to realize 
that the two equations ( |2.10|J2.lT| ) are identical to each other up to terms which have been 
neglected when the transport equation Q2.2|) was derived from the underlying exact equation. 
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Therefore, one might conclude that one can use one or the other equation to evolve the test 
particle coordinates and thus / in time. The difference between these two evolutions is in 
an order which was neglected anyway when deriving the transport equation from the exact 
one. Therefore, we now choose e.g. ( 2.11|) to evolve /. One has to realize, however, that now 



it is no longer the transport equation ( |2.2| ) one is solving. Instead, one can reconstruct from 
the test particle equation of motion that the corresponding transport equation is given by 



o = E 

i 

dp 



F (t) 

-Pi(t) 



- K(t,Pi(t)) 

F(t) 



1 - K(t,p) 



which clearly differs from (|2.2| ). However, by construction this new transport equation (|2.14|) 
is as close to the underlying exact equation as the original transport equation (|2T 



To summarize we have presented two ways to deal with the approximate transport equa- 
tion and the fact that the particle number is conserved by the exact equation. 

1. Extract a quantity which is exactly conserved by the approximate equation (here N). 
Make a test particle ansatz for the corresponding density (here /). Once the transport 
equation is solved one can also calculate / and N and figure out to which extent the 
conservation of N is violated in time. This provides a check for the accuracy of the 
approximation scheme which has served to derive the transport equation from the 
exact one. 

2. Make a test particle ansatz for the quantity which is exactly conserved by the exact 
equation but not by the transport equation (here /). Ignore one of the obtained 
equations of motion for the test particle coordinates. Usually one is tempted to ignore 
the equation obtained from the coefficient of the pure 5-function. Reconstruct a new 
transport equation from the test particle equations of motion. 

At first sight the second approach seems to be more appealing since there the recon- 
structed transport equation shares with the exact equation the full conservation of the par- 
ticle number. Unfortunately things are not always as simple as the chosen example seems 
to indicate. The problem is that in general it is not always possible to reconstruct such a 
transport equation from the test particle equations. As long as only the propagation of test 
particles is concerned (Vlasov-type equation) one can reconstruct a corresponding transport 
equation as shown above. If, however, also collisions of test particles are taken into account 
(Boltzmann-type equations) new problems arise. To see this we extend ( |2.2j ) by adding a 
collision term on the r.h.s. 

| [(1 - K (t,p))f(t,p)] + F(t) ^-f(t,p) = I coll [f] (2.15) 

with the additional condition 
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Jdplcou[f]=0. 



(2.16) 



We stress again that this equation has to be regarded as an approximate (transport) equa- 
tion. For the underlying exact equation we assume that it is too complicated to solve and 
that it conserves the particle number N. Again we observe that N is conserved by (|2.15 ) 



while N is not conserved. The most pragmatic (albeit doubtful) approach would be the 
following: One makes a test particle ansatz for /. For the propagation of the test particles 
between collisions one obtains again the overdetermined system of equations ( p,10| , |2.1lD . 
One disregards the first equation and uses the second one to evolve the test particle coor- 
dinates in time. One allows in addition for collisions between the particles according to the 
collision term on the r.h.s. of ( |2.15| ). The obstacle in that pragmatic approach is the fact 



that there is no corresponding transport equation at all which can be derived from ( [2.15 ). 
Reconstructing a transport equation from the test particle equation of motion yields 

Thus the collision term now is given by J co ii/ (1 — k) instead of I co \\. This new collision term, 
however, does not conserve the particle number 



Therefore, the reconstructed transport equation ( |2.17| ) has no advantages as compared to 



the original one ( [2.15 ). Using instead the original collision term J co n as described above is 



only a guess. It is not clear whether such an ad hoc postulated equation is as close to the 
exact equation as the transport equation (|2.15|) . 



Aiming at a field theoretical foundation for the description of transport processes we 
think it is not acceptable to have only a recipe for the test particle evolution which is not 
derivable from the underlying theory. If an already approximate transport equation is only 
approximately solved one loses the contact to the original exact equation more and more 
in a not controllable way. Therefore, we suggest to use the first approach described above, 
i.e. the test particle ansatz for /, since it solves the transport equation exactly. In addition, 
it provides the possibility to judge the approximation scheme which has led to the transport 
equation simply by evaluating the time dependence of the quantity which is conserved by 
the exact equation. 

After these precursory considerations concerning the use of the test particle ansatz we 
now turn to the derivation of a transport equation from the underlying field-theoretical 
Kadanoff-Baym equations. 



III. KADANOFF-BAYM EQUATIONS AND GRADIENT EXPANSION 

Following Kadanoff and Baym we start with the exact nonrelativistic equations of 
motion for the two-point functions 
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.d_ 



V 

2m 



D<(1,1') 



dl 
dl 



S ret (l, 1) D<(1, 1') + E<(1, 1) D av (l, 1') 
E ret (l, 1) D>(1, 1') + £>(!, 1) D av (l, 1') 



(3.1 
(3.2 



(3.3 
(3.4 



(3.5 
(3.6 

In ( |3.3| ) and throughout this work the upper (lower) sign refers to bosons (fermions) except 
where otherwise stated. The self-energies are connected via 



where we have introduced the two-point functions without ordering 

iD < (x,y) = ±(^(y)^(x)), 
iD > (x,y) = (4>(x)ip ] {y)) 

and the retarded and advanced quantities 

D* et (x,y) = 6(x - y ) [D>(x,y) - D<(x,y)} , 
D av (x,y) = e(y -x ) [D<(x,y) - D>{x,y)} . 



S rct ( 



x,y) ■ 



S HF (x,y) + 0(x o - y ) [E>(x,y) - E<(x,y)\ , 
: S HF (a;, y) + S(y - x ) [£<0r, y) - E>(x, y)} . 



(3.7 
(3.8 



The collisional self-energies S< and S> will be specified below. We have also introduced the 
time-local Hartree-Fock self-energy E HF (x, y) ~ 5(xo — yo)- In ( |3.1| , |3.2j ) the numbers denote 
short-hand notations for the space-time coordinates. 

In the following we will concentrate on systems where the dependence of an arbitrary 
two-point function or self-energy F(x,y) on its center-of-mass variable (x + y)/2 is weak 
(albeit not negligible). Therefore it is appropriate to introduce the Fourier transform with 
respect to the more rapidly oscillating difference variable x — y (Wigner transformation): 



F(X, p)= JdSi e ipu F(X + u/2, X - u/2) 
where F denotes an arbitrary two-point function or self-energy, 

F = D< £>>...,£<,.... 



(3.9) 



(3.10) 



From the definitions ( |3 . 3| . |3 . 4| . |3~5| . |3 . 6| . |3 . 9| ) we obtain the following transformation proper- 
ties with respect to complex conjugation: 



D<(X,p) 
D>(X,p) 
D TCt (X,p) 



-D<(X,p), 
-D>(X,p), 
D™(X,p). 



(3.11) 



Corresponding equations hold also for the self-energies. It is useful to introduce the real- 
valued quantities 



S<(X,p) = ±iD<{X,p) , S>(X,p) = iD>{X,p) 



(3.12) 
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and 



A(X,p) = -2lmD Tet (X,p) = 2lmD av (X,p) 

= i[& c \X )V ) - D^(X,p)} = i[D>(X,p) - D<(X,p)} 

= S > (X,p) T S < (X,p). (3.13) 

The quantity A is the spectral function. On account of the (anti-)commutation relation for 
bosonic (fermionic) fields 

^(t, x) ^(t, y) =F ip\t, y) x) = 5(x - y) (3.14) 

the spectral function is normalized to one: 

J^A(X,p) = l. (3.15) 

The quantity S K can be interpreted as a generalized phase-space density. The particle 
number can be obtained from S K as 



N(t) = Jd A x j--^S < (t,x;p). (3.16) 

For a system which has reached thermal equilibrium S < (X, p) does no longer depend on X 
and is given by 

S^(p)=n B , F (p )A(p) (3.17) 

where Ub,f is the thermal Bose (Fermi) distribution. 

For a general off-equilibrium situation we are aiming at equations of motion for the real- 
valued quantities S K , A, and Re/) rct = ReD av . (Note that S l> can be obtained from 5' < 
and A according to ( |3.13|) ). In principle, KeD iet can be obtained from A via a dispersion 
relation which can be easily obtained by Wigner transformation of the definitions (|3.5|,p.6|) : 



ReD^(X,p) = -°V -A(X; k ,p) (3.18) 

J 2tc p — Kq 

where V denotes the principal value. It is easier, however, to derive the equation of motion 
for ReZ) ret from an equation of motion for D ret . Using ( |3.1| , |3.2|) and the definition (|3.5| ) we 

get 

(i-^ + ^) ^ rct (l, 1') = 5 (4) (1, 1') + Jdl S rct (l, 1) D re \l, 1') . (3.19) 



As we will see below the final result for ReD ret agrees with (|3.18j ). 

As already mentioned we are interested in systems where the dependence of the two- 
point functions on X is weak. This allows to Wigner transform the equations of motion 



( |3.1| , |3.2| , |3.19| ) and to neglect all terms with more than one derivative with respect to the 



center-of-mass variable X. We get the following equations of motion in first order gradient 
expansion: 
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-n 
P 



Po — ~~ ] D<(X,p) - - [p - D<(X,p)] 



2m 



2m 



£ rct (X,p) L><(X,p) + E<(X,p) L> av (X,p) 
- ~ [£ ret (X,p), 5<(X,p)] - % - \E<{X,p), D™(X,p)\ 



(3.20) 



/r 2 



po - f - ] 5>(X,p) - 1 [p - £»(X,p)] 



2m 



2m 



£ rct (X,p) D>(X,p) + E>(X,p) L> av (X,p) 
- I [£ ret (X,p), D>(X,p)) - I [£>(X,p), L> av (X,p)] 

(po-^)^(X,p)-l[p -|-,^(X,p)]^ 

i + s rct (x,p) ^ rct (^p) - ^ [s rct (^,p), ^ rct (^p)] 

where we have introduced the generalized Poisson bracket 

[A, B] = d Xo Ad po B - d Po Ad Xo B - V X A V P B + V P A V X B . 



(3.21) 



(3.22) 



(3.23) 



It is worth noting that the drift terms typical for transport equations show up in the ex- 
pressions 



-2 -* 

[po - |-, B] = -d Xo B -^-V X B 
2m m 



(3.24) 



We note that for the derivation of ( |3.20| - |3~2"2"D we have not assumed anything about the 
derivatives with respect to p. This is different to the approach of [0] where systems were 
studied for which mixed double derivatives like V ' x V p A are also negligibly small. We will 
come back to that subtle difference below. 

Next we extract from ( p. 20 - 3. 22 ) real and imaginary parts using ( p. 12 , 3. 13 ) and the 
decompositions 



£ ret/av (X,p) = Re£ ret/av (X,p) =p ^iT(X,p) 



(3.25) 



where we have introduced the (real-valued, non-negative) width 

T = z(S> - £<) = z(S rct - £ av ) . (3.26) 
In ( |3.25| ) the minus (plus) sign refers to the retarded (advanced) self-energy. We find 

po - i J S< = ReS rot S< ± iS< ReD TCt - i[T, S<] + ^[±^<, A] , (3.27) 



(po - 


2m , 


) S< 


[Po - 


p A 
2m' 


S<] 



Sretl 



(3.28) 
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[Po - — I A = ReS rot A + T R eJ D ret , (3.29) 
\ 2m y 

[p - = [ReS ret , .A] + [r, ReD ret ] , (3.30) 

( Po ~ 2m") R6jDrCt = 1 + Re ^ r6t ReDCt 'l^^-^^' ^ 3 - 31 ^ 

- 1 

[p - |— , ReZ) rot ] = [ReE ret , ReD ret ] - -[T, A] . (3.32) 



Note that «S < is real- valued (cf. Eq. (ft.ll ) and the remark concerning self-energies). 



Obviously using the gradient expansion up to first order we have obtained six equations 
for three quantities (S < , A, and ReZ) ret ). Therefore some of these six equations are either 
redundant or provide constraints which have to be fulfilled by the dynamical quantities to 
make sure that their evolution is in accordance with the first order gradient expansion. To 
say it in other words: Even if one starts out with a configuration (represented by the two- 
point functions at initial time) for which the dependence on X is weak, a peculiar choice for 
the self-energies might drive the system out of the regime where the gradient expansion is 
valid. Thus it would not be surprising if the gradient expanded equations of motion provided 
constraints on the dynamical quantities. 

As we will see now there are indeed redundant equations as well as constraints. Let 
us start with the purely algebraic equations ( |3.29| , |3.3lD . They already yield the following 



expressions for two of the three quantities of interest: 

A(X, p) = _ r(X,P) - 2 , (3.33) 

(po-£-Re£-t(X,p)) +\Y\X lP ) 

ReD-(X, P ) = *-£-B^(*,P) 

(po-fi-ReE-t(x,p)) +\T%X,p) 

For the retarded and advanced two-point functions this yields 

D Ict/av (X,p) = =3 = . (3.35) 

Po-£-£ rct / av (^P) 

These relations are very well known for systems which are completely homogeneous in space 
and time, i.e. where the zeroth order in the gradient expansion is sufficient. It is interesting 
to observe now that these relations are still valid in first order gradient expansion |15|,|TT . 
We also note that inserting ( |3.33|) in ( [3.18| ) yields (|3.34j ) provided that the used self-energy 
is indeed retarded, i.e. satisfies the dispersion relation 

-r(X; (3.36) 

Finally one can show that the expression for the spectral function ( |3.33| ) indeed fulfills the 
normalization condition (|3.15|) . This is discussed in Appendix 
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Next, we observe that the expressions ( |3.33| , |3.34D identically solve (|3.30| , |3.32|) . To see 
this one might simply insert ( p. 33 |T34 ) and perform the somewhat tedious brute-force cal- 
culation. Instead we will take here the following elegant way |I9|| : We combine ( p.30| , |3~32| ) 
to a complex equation. ( |3.32|) — ^( p.30|) yields 



D 



rctl 



0. 



(3.37) 



It is trivial to see that the expression ( |3.35| ) solves this equation. It is interesting to note 
here that the same reasoning holds for relativistic scalar bosons. For relativistic fermions, 
however, things might be more complicated since in this case the generalized Poisson bracket 
of the matrix valued quantities D ret and (D^)" 1 may not vanish any more [19]. To the best 
of our knowledge this subtlety has never been worked out in detail (concerning a closely 
related case cf. p0|). 

We conclude that the equations ( |3.30|J3.3^ ) provide no additional information as com- 
pared to ( |3.29| , |3.3lD ; they are redundant. Therefore, the set of equations ( |3.27| - |3.32"D contains 
at most four independent equations for the three quantities of interest. 

Now we turn our attention to ( |3.27| , pT2"SD . Before analyzing them for the general case 
of arbitrary width we briefly discuss their meaning for the quasi-particle regime: In this 
limit the gradient terms in ( |3.27|) and also the second term on the r.h.s. vanish and this 
equation becomes a purely algebraic (constraint) equation demanding that S K ought to be 
proportional to an on-shell ^-function. In other words this equation reduces to the mass-shell 
constraint discussed e.g. in 



The traditional on-shell transport equation is recovered 
from ( |3.28|) by integrating over the energy ||. As already pointed out the l.h.s. provides the 
drift terms. The first two terms on the r.h.s. yield the collision terms as will be discussed 
in more detail below. The third term on the r.h.s. gives the Vlasov contribution. The last 
term on the r.h.s. vanishes in the quasi-particle limit. Up to now the role of this term in 
an off-shell transport theory is not clear. However, once the equations of motion for test 
particles will have been written down the meaning of this term will become obvious. 

Let us now come back to the more general case of arbitrary width. Both equations 
( |3.27| , j3!2~8D mix zeroth and first order gradient terms. We first concentrate on the zeroth 
order terms: 



-o 

p 



Po- — - Re£ ret S< = ±iZ < ReD ret + o(d x ) 
2m J 

TS < = ±*S< A + o(d x ) ■ 



Since we already know expressions for ReZ) ret and A we insert ( |3.33| , |3.34D and obtain 

po - £ - Re£ ret 



po- — - ReS ret ) S< = ±zS< 
v 2m I 



Po 



£-ReS-) 2 + IP 



Re£ r 

r 



Po-£-ReS-) 2 + in 



+ o(d x ), 
+ o(d x ). 



(3.38) 
(3.39) 

(3.40) 
(3.41) 



At first glance these two equations seem to contain the same information in zeroth order 
gradient expansion: 
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* =" ±tS _ N2 . »(<M (3. 1.2) 



obtained from (|3.40|) by dividing by (Po ~ ~ Re£ ret ) and from (|3.41|) by dividing by T. 
However, we have implicitly assumed in these steps that both factors are not too small, i.e. 



2m 

1 



Po-h- Re£«* 



r 



o(d x ) = o{d x ) , (3.43) 
o(d x ) = o(d x ) • (3.44) 



This can only be true if we are not in the quasi-particle regime. There, the width is very 
small. Therefore, e.g. (|3.44j ) might be invalid. Fortunately we are here interested in the case 
where the width is finite.Q Thus, ( 3.42|) holds and provides information about the size of a 



specific combination of Green functions and self-energies: 

±iS< -dS < = o(d x ) (3.45) 

with 



d=(p Q -^-- ReL mt \ + ^r 2 . (3.46) 



Note that d is the denominator of A as well as ReZ) ret . Since all our consideration are 
based on the assumption that space-time derivatives are small we conclude that the l.h.s. of 
( |3.45| ) also has to be small. This is a consistency condition which might be checked in actual 
calculations. Since we keep terms linear in the gradients the l.h.s. of ( |3.45 ) in general cannot 



be neglected. However, when it appears in gradients we can neglect such a combination, 
i.e. for arbitrary B: 

[±iE< B] = [dS K , B] + [±«S< -dS < , B] = [dS<, B] (3.47) 



where we have neglected the gradient which is effectively of second order in the gradient 
expansion 0. Therefore, it turns out to be consistent with the gradient expansion up to 
(including) first order to replace izZ^ in the first order gradients in ( |3.27| , |3.28D by dS < . 
Then, after some straightforward manipulations we find that both equations reduce to the 
same transport equation (cf. [[K]]) 

\TA [po - - ReS ret , S<] - l -A [r, (Po - ^ - ReS ret )S<] = Y S< T i^A . (3.48) 



Nonetheless, in the end we want to have an equation which is also valid for small width, i.e. shows 
the correct quasi-particle limit. Indeed, as we shall see below, this requirement is fulfilled. Since 
a transport equation for the quasi-particle limit is well-known we can treat in our derivation the 
case where the width is large. 
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With this equation we have reached our goal to derive in first order gradient expansion 
consistent equations for the three quantities S K , A, and ReD ret . These equations are given 
by ( |3.48| ), ( |3.33| ), and ( |3.34j ), respectively. They have also been derived in [fT6|JT7f1 . It is 



interesting to note that only the equation for S < is a differential equation. 

Let us recall again that the derivation of this latter equation has been obtained in a 
two-step procedure. The first step was the formal gradient expansion yielding two equations 
for 5 ,< , namely ( |3.27| ) and ( |3.28| ). The second step was to realize that a special combina- 
tion ( |3.45| ) of self-energies and Green functions is effectively of first order in the gradient 
expansion. Without that second step the two equations ( |3.27| , |3.28"D would not yield identical 
results and it would be unclear which equation one should use to evolve S K in time. Of 
course, it is interesting to find out for which physical situations (|3.45|) is realized since in 
turn this indicates where the gradient expansion is justified. For that purpose we rewrite 
condition (|3.45|) in the way as it appears in the traditional transport equation (|3.28|) : 



rS<TiZ<A = rA [^--^1] =o{d x ). (3.49) 




Obviously, there are two cases when this special combination of self-energies and two-point 
functions is small. The first case is the traditional quasi-particle regime. There the width V 
is small and from 

r = o(d x ) (3.50) 

one might justify the gradient expansion.^ The second case, relevant for our purpose of 
describing states with potentially large width, concerns situations where the difference of 
the fractions in the brackets in ( |3.49| ) is small. Actually these two fractions define off-shell 



generalizations of phase-space densities, one given by the two-point functions 

A{X,p) 

and one demanded by the self-energies 

MX,P) = ^^. (3.52) 

It is the difference of these two functions which in general determines the time evolution of 
the system under consideration (see also e.g. [15,11,161). The smallness of this difference is 



an alternative possible justification of the gradient expansion if the width of the states is not 
small. Indeed, this difference is small near thermal equilibrium since there the collisional 
self-energy is given by [0 



2 Naively one might suspect that the spectral function becoming sharply peaked invalidates the 
given argument in the region of interest, i.e. near the on-shell point. However, in this case one 
should rather compare energy integrated expressions. Then upon integration over the energy the 
spectral function in front of the brackets in ( |3.49| ) is basically replaced by one. 
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± iS < (p) = n B ,F(po) r(p) 



(3.53) 



while S< is given in (|3.17 ). This shows that both n and n^, approach the same thermal 
distribution. The fact that ( |3.49| ) is fulfilled near thermal equilibrium does not automatically 
imply that it is not fulfilled outside of this regime, i.e. far away from thermal equilibrium. 
Whether ( |3.49| ) is realized or not for a given physical situation has to be checked in actual 
calculations. If it was violated the consequence would be that the process at hand cannot 
be described by first-order gradient expanded equations. 

The transport equation ( |3.48| ) is the main result of this section. To gain more insight 
into this equation we first rewrite the r.h.s. in a more common form. Using ( |3.13|J3.2B| ) we 
get 



YS < t iX<A = =F iZ < S > 



(3.54) 



which is the well-known collision term, the difference between loss and gain terms. In the 
test particle description we are aiming at, the l.h.s. of (|3.48 ) describes the propagation of 
the test particles. For a first interpretation of these terms on the l.h.s. we study the quasi- 
particle limit. In spite of the fact that the derivation of ( |3.48| ) is doubtful for infinitesimally 
small width (as outlined above) we will find that ( |3.48| ) nonetheless reproduces the correct 
transport equation for the quasi-particle regime. In this regime one may define an on-shell 
phase-space density n via 



S<(X,p) = n(X;p = e(X,p),p)A(X,p) 
where e is the on-shell energy, i.e. the solution p = e of the equation 



Po 



-e> 

2m 



-ReS rct (X,p) = 



(3.55) 



(3.56) 



In this way ( p.48|) can be transformed into a transport equation for n. We integrate the 
transport equation over po and note that for vanishing width both A and ^IVl 2 approach 
the same 5-function: 



A 1 



P 2 



2n 5(p - — - Re£ ret ) = 2vr z 5(p - e) for T 
2m 







(3.57) 



with the renormalization factor 

z(X,p)= ( 



1 - 



<9Re£ ret (X,;p)' 
dp 



(3.58) 



P0=e 



The collision term, i.e. the r.h.s. of ( |3.48| ), becomes after integrating over p Q : 

- z (±iE<(l±n) - if^rij 
The first term on the l.h.s. of ( |3.48| ) yields the Vlasov part (see also ( |3.24| )): 



(3.59) 



PO=e 
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[Po - e, n] 



(3.60) 



PO=e 



Note that the ^-factor is implicitly incorporated here on account of 



ret 



(3.61) 



P0=e 



and corresponding equations for derivatives with respect to momenta. The second term on 
the l.h.s. of (|3.48|) is a genuine off-shell contribution. After performing the integration over 
Po this term vanishes for r — > 0. The origin of this term can be traced back to the last term 
on the r.h.s. of ( |3.28| ). As already emphasized the role of this term in an off-shell transport 
theory still requires a proper interpretation which we will give below in the test particle 
picture. To summarize, we have recovered from ( |3.48| ) the traditional Boltzmann equation 
in the quasi-particle limit. 

In the next sections we will elaborate on the more general case of non- vanishing width. 
Among other things we will address two questions: 1. Is the particle number still a conserved 
quantity when the transport equation ( |3.48 ) is used to describe the time evolution of S' < ? 
2. What is the role of the off-shell contribution in (|3.48|) ? 



IV. EFFECTIVE PARTICLE NUMBER 



We are aiming at a test particle representation for the transport process described by 
(p. 48]) . For a test particle representation to make sense it is mandatory that the represented 
quantity is conserved as outlined in Sec. [TI[ Suppose that we have chosen the self-energies 
such (conserving approximation P,|7|,|T6ll ) that the particle number ( |3.16|) is conserved for 



the full quantum theory defined by the exact equations of motion ( |3.iy3.2j ). This does not 
necessarily mean that also the transport equation (|3.48| ) conserves this quantity. Indeed, as 
we will argue below this is not the case. The reason is simply that by the coarse graining 
process carried out by the gradient expansion the part of S < (X, p) which is highly oscillating 
in X has been neglected. Of course, the temporal deviation in the particle number is 
effectively of higher order in the gradient expansion since up to first order the exact and the 
approximate equations of motion carry the same information. However, it would be very 
unpleasant if one had to represent a quantity by test particles which is not fully conserved 
on the level of approximations. Thus, after showing that the particle number is not fully 
conserved we shall search for another quantity which is exactly conserved by the transport 
equation. 

To study the time evolution of the particle number ( |3.16| ) we divide ( |3.48| ) by \TA and 



integrate over the space coordinates and the four-momentum. All gradients with respect to 
x and p vanish by partial integration (and neglecting surface terms). Finally we get 



±N(t) = 4- [d 3 x f-^-S<(t,x;p) 
dt y ' dtJ J (2tt) 4 v ' ,yj 



d fjL. f #P c <^ , fjL. f d4 P 2 c ;v c 
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with 



= dR ^(X,p) + p -^-Re^\X,p)dT(X,p) ^ 

Note that the quantity 1 — K is the off-shell generalization of the inverse quasi-particle 
renormalization factor 1/z as defined in (|3.58 ). 



Now we have to specify the self-energies. To keep things simple we choose a two-body 
interaction 

^ dV V f (t, x) ^(t, £') - x*\) ip(t, x') ip(t, x) (4.3) 
and evaluate the collisional self-energies in the Born approximation 

?S (X ' p)= i(2^(2^(2^ (27r) ^(P + Pi-l*-**) 

x^(^(p-P2)±^(p-P3)) 2 5 < (X,p 1 )5 > (X ) p 2 )5 > (X,p3), (4.4) 

±ZS (X ' P) = 7 (2^(2^(2^ (2?r) 5 + -P2 -Ps) 

x^(^(p-P2)±^(p-P3)) 2 5 > (X,p 1 )5 < (X,p 2 )S < (X,p 3 ). (4.5) 



For a more elaborate treatment of self-energies we refer to 

The two contributions on the r.h.s. of ( |4.1|) show completely different structures. Es- 
pecially we note that the first term belongs to the part of the transport equation which 
describes the propagation of modes while the second term comes from the collisional part. 
It is hard to conceive that these two terms can cancel each other. Indeed, we have not 
managed to rearrange these terms in a way that one can get an idea why the r.h.s. of ( fi.lp 
should vanish. We therefore strongly conjecture that N(t) is not conserved by the transport 



equation ( 3.4£| ). At this stage a comparison to the approach of []7| is in order. There terms 



like the first one on the r.h.s. of ( |4.1| ) are neglected since it contains a double derivative 
- one with respect to t and one with respect to po. We cannot neglect such contributions 
since we have only assumed that space-time derivatives are small. We have not made any 
assumptions about the smallness of derivatives with respect to momenta. Even more impor- 
tant, we are not aiming at a quantity which is nearly conserved by the transport equation 
up to terms which are of higher order in the gradients. To represent the transport equation 
by test particles we need a quantity which is exactly conserved. Of course, this difference 
is somewhat subtle since e.g. the particle number is exactly conserved by the full quantum 
mechanical equation of motion. However, it is not the latter equation we want to solve by 
a test particle ansatz but the transport equation obtained by gradient expansion. 

We will now construct a quantity which is conserved by (|3.48|) . The first step is to get 
rid of the collisional part. By inspecting (f4.4|^4.5| ) we observe that the following integral 
vanishes 

(±zS < 5 > -zS>S<) = 0. (4.6) 



(2*)' 



16 



This suggests to directly integrate (|3.48|) over momenta and space coordinates instead of 
dividing it by ^TA first. Indeed, for the quantity 



S = -TAS < 
2 

we derive from ( |3.48D the following equation of motion (see Appendix |B|): 

[po - |" - ReE ret , S] - - [r, (p - f - - Re£ ret )<S] = z£>S< =f *S<S> 
2m 1 2m 



Thus we find 



(2vr) 4 dt 



W) 1 



SK. 



(4.7) 



(4.8) 



(4.9) 



We conclude that the space and momentum integral of the quantity 

1. 



S< = S (1 - K) = 2 r ^ 5< i 1 - K ) 



(4.10) 

is conserved by the transport equation ( |3.48| ): 

-S<(t,x;p) = Jd 3 x J-^-(±iE<S>-iE>S<)=0. (4.11) 



(27T) 



In the following we will refer to N as the effective particle number. In the quasi-particle 
limit, r — > 0, we find 



S < (X,p)^n(X;e,p)2nS(p -e) 
where we have used ( |3.55|J3.57| ) and 

1 - K z' 1 . 

Thus, in the quasi-particle limit there remains an interesting difference between 

S < (X 1 p)^n(X;e,p)27rz5(p -e) 



(4.12) 



(4.13) 



(4.14) 



and S < . The quantity which "counts" the quasi-particles is indeed S < while the renormal- 
ization factor z appears in S < which takes into account that there is some strength in the 
spectral function (far) away from the quasi-particle pole. 

As a second limiting case we evaluate S < for thermal equilibrium. Using (3.17) we get 



3£(p) = "b,f(pd) \T{p)A\p) (1 - K{p)) . 



(4.15) 



For constant width and vanishing ReS ret the quantities and are depicted in Fig. [I]. 
Obviously the function |r^4 2 is more strongly peaked than A (also cf. |12|). The asymmetric 
forms are caused by the Bose enhancement. 
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Comparison of ( |4.15| ) with ( |3.17|) suggests to introduce — also for non-equilibrium situ- 



ations — an effective spectral function 

A{X,p) = h(X,p)A 2 (X,p) (1 - K{X,p)) . (4.16) 

Indeed, as will be shown in Appendix [A] the effective spectral function is also normalized 
(cf. (|3.15|) ). Thus we have a correspondence of the quantities S < , N, and A of the full 



theory with the quantities S < , N, and A of the coarse grained theory obtained by gradient 
expansion. However, a word of caution is in order here. For the interpretation of A as an 
effective spectral function (and also of N as an effective particle number!) it is mandatory 
that A is always non-negative. This may not be the case for arbitrary self-energies. Since 
T is always positive one has to make sure in actual calculations that 1 — K is non-negative. 
Otherwise a test particle ansatz would not make any sense. Finally we note that A was also 
introduced in |TjJ as the spectral information which enters the entropy density. 



V. EQUATIONS OF MOTION FOR TEST PARTICLES 

The next step is the test particle ansatz for S K and the derivation of the equations of 
motion for the test particles from the transport equation ( |3.48j ): 

S<{t, x- p) ~ * (3) (* ~ 5 (Po - Ei(t)) 5 {3 \p- Pi{t)) • (5-1) 

Note that the energy Ei of the test particle i is a free coordinate, not restricted by a mass-shell 
condition. Instead of ( |3.48| ) we use the equivalent equation ( ^4.8| ). While the r.h.s. describes 
the collisions the l.h.s. will give the propagation of the test particles between collisions. To 
get the equation of motion for a single test particle between two collisions we have to insert 

(5.2) 

into 

S] - - [r, (po-~ ReS ret )5] = . (5.3) 
In general this yields an equation of the following type 

Y (a(t,Xi; Ei,pi) + b(t,Xi] Ei,pi)d Po + c(t,Xi] E i} pi)V x + d(t,Xi] E h p i )V i 

i 

®5^{x-x i )8{p Q -E i )8®{p-p i ) = U. (5.4) 

A solution for this equation is obtained by demanding that all coefficient functions a, b, 
c, and d have to vanish. In general this gives eight equations for the seven test particle 
coordinates E i: x i} and f>i. If and only if the quantity (in our case S ,< ) which is represented 
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by test particles corresponds to a conserved quantity (here N) then the coefficient a in (|5.4|) 
vanishes and one only has to fulfill seven equations instead of eight. This is the reason 
why we have insisted to find an exactly conserved quantity. Otherwise one would have to 
deal with an overdetermined system of equations. Neglecting one of the obtained equations 
would not yield a solution for the transport equation one actually wants to solve (cf. the 
corresponding discussion in Sec. |J). F° r t ne case at hand we get the following equations of 
motion for the test particles 

TP. _ P? "D^Vret \ 

(5.5) 





^ + V Pi Re£ ret + ^ % — V P T I = c(t, x % - E^fr) , (5.6) 



jjj. _ £i ReE ret \ 

V^ReS ret + 1 2m T V X T =d(t,x l ;E l} p i ). (5.7) 

Newton's equations of motion are obtained by neglecting K, i.e. putting the renormalization 
factor 1/(1 — K) to 1, and disregarding all terms where the width T enters. In this limiting 
case one observes that the energy of a test particle is only changed if the real part of 
the self-energy, i.e. the classical potential, depends explicitly on time. If the potential 
shows an explicit momentum dependence the classical expression for the velocity is modified 
accordingly. For finite width but neglecting the energy and momentum dependence of the 
self-energies ( j5~5"| - |5.7| ) basically reduce to the test particle equations presented in []T7|] . The 
only remaining difference is that there a relativistic system is studied. 

Let us now come back to the full equations of motion ( |5.5| - |5.7D . Obviously the terms 
which involve the width yield genuine off-shell contributions to all equations of motion since 
on the mass-shell the combination Ei — J-*- — Re£ ret vanishes by definition (cf. fl3.56| )). It is 
interesting to calculate the time evolution of the off-shellness of a test particle defined by 

AEi(t, x t ; Ei,pi) =Ei-^-- Re£ ret (t, E u p % ) . (5.8) 



From the equations of motion (|5.5| - j5"77 ) we find 



This equation has also been presented in [17] and in [18|. To the best of our knowledge 



the full set of equations of motion ( |5.5| - |5.7D has never been derived before. Obviously the 
time evolution of the off-shellness is caused by the respective last term in the brackets on 
the l.h.s. of ( |5.5| - |5.7|) . These contributions can be traced back to the second term on the 



l.h.s. of ( p.48| ). This clarifies the meaning of the latter term for an off-shell transport theory: 
It provides the time evolution of the off-shellness. Actually this off-shell contribution in the 
transport equation ( |3.48j ) is caused by the respective last term on the r.h.s. of ( |3.27| ) and 



( |3.28j ) via the replacement (|3.47 ). Without that latter replacement the meaning of these 



terms in ( |3.27| , |3.28"D would have been completely unclear. 
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The equations of motion ( |5.5| - |5.7| ) obeyed by the test particles between collisions are the 
main result of this section. It is useful to show that they reproduce the correct quasi-particle 
limit. The quasi-particle energy of a test particle is defined as the solution = e^t^x^Pi) 
of the equation 



-f) 
2m 



Re£ ret =0. 



Differentiating (|5.10| ) e.g. with respect to Pi yields (also cf. ( p.61|) ) 



(5.10) 



v r 

Vl 1 1 - K(t,Xi] e i: pi) \m 
Thus we find the expected relations 



^ + V Pi ReZ Tet (t,x i ;e i ,p 



X i 

Pi = -V^ei . 

The time evolution of the on-shell energy is governed by 
1 



(5.11) 



(5.12) 
(5.13) 



1 - K(t,Xi, €i,Pi] 



d t ReY? et {t } Xi]ei,pi) = z(t, x { \ e h pi) <9 t ReS ret (t, x { ; £i,Pi) . (5.14) 



Finally we comment on the evaluation of the self-energies which enter the equations of 
motion ( |5.5| - [5.7| ) and of course also the collisional part on the r.h.s. of the transport equation 
( |3.48| ). Except for the Hartree-Fock self-energy which is comparatively easy to incorporate in 
transport calculations all self-energy contributions can be traced back to the determination 
of ±iE < and iE > via the relations ( |3.26|j3~l36|j4.2f ). In the Born approximation to two-body 
collisions these self-energies are given in ( |4.4| , [4.5| ). These quantities should be expressed in 
terms of S < since the latter quantity is represented by test particles according to (|5.1|) . We 
find for the Born rates 



iT>{X,p) 



(2tt) b { >{p + pi -P2-P3) - {v{p-p2) ±v{p-p 3 )) 



(2tt) 4 (2tt) 4 (2tt) 4 



A 2 ± S< 



A* ± S< 



lT 1 A 1 {l-K 1 ) \T 2 A 2 (l-K 2 ) ir 3 ^ 3 (l-^3) 



(5.15) 



±iZ<(X,p) 



° 1 1>2 " ' ''" (2n) 4 5 (4 \p + p 1 -p 2 -p 3 )]-(v(p-p 2 )±v(p-p 3 )) 2 



(2tt) 4 (2vr) 4 (2tt) 4 
A 1 ± S< 



S< 



5 3 < 



ir x Ai (1 - K x ) \Y 2 A 2 (1 - K 2 ) ir 3 ^3 (i - K 3 ) 



(5.16) 



with = S < (X,pj) etc. Obviously the fact that S < instead of S < is represented by 
test particles causes rather non-trivial modifications for the collision integrals. In turn, 
disregarding these modifications amounts to solving transport equations which might be 
rather different from the ones one actually wants to solve. 
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VI. SUMMARY AND OUTLOOK 



In this work we have presented a derivation of a transport equation making use of the 
gradient expansion but without utilizing the commonly used quasi-particle approximation. 
Therefore, the derived equations can be used for the description of processes where the width 
of the involved modes might be arbitrarily large.0 In principle, the gradient expansion yields 
six equations for the three quantities of interest. We have successively shown which of these 
equations are redundant by various rearrangements. In the course of this procedure it has 
become apparent that a specific combination of self-energies and Green functions — which 
is seemingly of zeroth order in the gradient expansion — must be treated as being of first 
order. This identification ( p.45|) is mandatory to derive a contribution to the transport 
equation which finally leads to the description of the time evolution of the off-shellness of 
the test particles. As a next step a quantity — the effective particle number — has been 
isolated which is exactly conserved by the transport equation. This has opened the way to 
a test particle representation of the density S K which corresponds to the effective particle 
number. The equations of motion for the test particles have been presented. Finally the 
Born collision rates are evaluated within the test particle representation. 

Remarkably the effective particle number is not identical to the particle number which 
is exactly conserved by the fully quantum field theoretical equations of motion. Of course, 
the particle number is approximately conserved by the transport equation since up to first 
order in the gradient expansion the approximate (transport) and the exact (quantum field 
theoretical) equation carry the same information. To apply a test particle ansatz to the 
transport equation one needs, however, a quantity which is exactly conserved. Nonetheless, 
in actual simulations the particle number N of the full theory might serve as a test for the 
accuracy of the gradient expansion. If during the time evolution N deviates drastically from 
its initial value the gradient expansion has to be regarded as inappropriate for the system 
one wants to simulate. 

Throughout this work we have strictly distinguished between the (possibly off-shell) 
propagation of the test particles between collisions on the one hand side and the binary 
collisions on the other. This distinction might seem to be artificial since the width which 
allows for an off-shell propagation is caused by the collisions. Since our equations are derived 
from the underlying quantum field theory in a well-defined and controlled way we do not have 
to be afraid of double counting. The physical picture we have in mind which corresponds 
to our approach are collisions which are so frequent that particles do not come back to 
their mass-shell before the next collision happens. Therefore part of the collision process 
has to be incorporated in the propagation by allowing the propagation of test particles with 
arbitrary energy. An additional aspect which we have not touched here but is in principle 
straightforward to include is the treatment of resonances with large decay width. Here it 
is very natural that a propagating resonance is simulated by a bunch of test particles with 
variable energies. We note that our formalism treats collisional and decay width on equal 



Of course, for a non-relativistic formalism to make sense one has to assume that kinetic and 
potential energy and also the width is small as compared to the rest mass of the particle. 
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footing. 

We would like to comment briefly on the possible numerical realization of the presented 
approach. To calculate for a given space-time point the spectral function and the width 
which enter the collision integrals ( ^.15| , ^T6|) one has to know the self-energies for arbitrary 



four-momenta. This requires the use of much more test particles as compared to simulations 
which are restricted to the quasi-particle regime. In present simulations of off-shell effects 
in transport theories rather ad-hoc recipes are in use to overcome this problem. In [I7| the 



self-energies are not determined independently for every point in space but rather from an 
average over a larger volume. Whether this is a valid approximation to the full treatment of 
the self-energies remains to be seen. In [|18j the self-energies for a given space-time point are 
parametrized by an expression adopted from the case of thermal equilibrium. The thermal 
parameters are locally chosen such that a reasonable fit on the shape suggested by the test 
particles is obtained. It has been checked there that time as well as momentum averages 
of the parametrized self-energies agree with the corresponding results from the test particle 
representation. 

Next we compare our formalism to recent other approaches. In it has been worked 
out in great detail how general self-energies have to be chosen such that particle number and 
energy are exactly conserved by the full quantum field theoretical Kadanoff-Baym equations. 
It has not been shown there that these quantities are exactly conserved by the approximate 
gradient expanded equations. As we have outlined above we conjecture that this is indeed not 
the case. To some extent it is trivial that quantities which are exactly conserved by the exact 
equations are (at least) approximately conserved by the approximate equations. Nonetheless 
it is interesting to see how this approximate conservation comes about in the approximate 
transport formalism J7j. In our approach, however, we had to answer a somewhat different 
question. For the test particle realization we had to find a quantity which is exactly conserved 
by the approximate transport equation. 

The purpose of ]T7|] and |TB[ was also to present a test particle description of off-shell 



transport processes. In spite of the fact that in |T7j the transport equation (|3.48 ) was de 



rived in basically the same way as presented here, the approach of |L7j is very different in 
spirit. There the gradient expanded Kadanoff-Baym equations were only used to determine 
the evolution of test particles between collisions. In contrast, the information about how 
collisions between these test particles have to be treated was taken from a completely dif- 
ferent source. Since in |l7j it was anyway not the purpose to solve equation ( |3.48D it did 



not matter whether the test particle equations of motion were derived from a test particle 
ansatz for S < or for S < . Obviously, in the approach presented here we take the transport 
equation ( |3.4ij| ) serious as being an equation which can be obtained from the underlying 
quantum field theoretical equations in a straightforward and well-controlled manner by gra- 
dient expansion. A second difference between our approach and the one presented in [17 



is the fact that the energy-momentum dependence of the self-energies was neglected for the 
test particle equations of motion. Therefore, the test particle equations presented here are 
more general than the ones given in ||17|| Finally we studied here a non-relativistic system 



4 After finishing this work a preprint [21| appeared where the authors of [O] generalized their 
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while in [fL7|] the case of relativistic bosons was treated. 

In [18] the starting point was the transport equation ( 3.28| ) but neglecting the last term 
on the r.h.s. The equation ( [5.9| ) which describes the time evolution of the off-shellness 
was introduced by hand based on physical plausibility arguments. The collision terms were 
treated in the way adopted from on-shell transport theories by simulating S K instead of S < 
by test particles. 

While discussing at length the (non-)conservation of the particle number we have not 
touched the issue of energy conservation. In view of the fact that it was non-trivial to 
find a conserved effective particle number it might be no surprise to realize that also the 
conservation of energy as defined for the full quantum field theoretical equations becomes a 
problem for the gradient expanded equations. To clarify the important question of energy 
conservation is beyond the scope of the present paper. We only want to note here that 
the energy as defined for the full theory (see e.g. [§,[7|,[l6j) has to be at least approximately 
conserved by the transport equation. In principle this provides a check on the validity of 
the gradient expansion. 

Finally we give a brief outlook on the generalization of our approach to relativistic 
systems. Most of the presented formalism can be immediately generalized, especially for 
scalar particles. As already mentioned above (after Eq. ( |3.37| )) it might appear for fermions 
that the set of four equations for spectral function and real part of the retarded propagator 
derived in first order gradient expansion cannot be trivially reduced to two independent 
equations. This might lead to new constraints for self-energies and Green functions. A 
second potential problem of a relativistic off-shell transport theory concerns the possibility 
of space-like modes. In principle, the spectral function might contain strength in the sector 
where the momentum is larger than the energy. Thus it can happen that in collisions space- 
like test particles are created which travel faster than the speed of light. This is clearly 
an unpleasant feature of a theory which should respect causality. For relativistic systems 
one therefore has to think about the proper treatment of the space-like part of the spectral 
function. 
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relativistic test particle equations by also including energy-momentum dependences of the self- 
energies. 
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APPENDIX A: NORMALIZATION OF SPECTRAL FUNCTION AND 
EFFECTIVE SPECTRAL FUNCTION 



In this appendix we will prove the normalization of the spectral function A as given in 
( p.33| ) and the effective spectral function A as defined in ( |4.16| ). We implicitly assume in 
the following that A is always non- negative (cf. the remarks after ( |4.16|) ). For the spectral 
function we find 

oo oo 

/^(X,p) = -2Im /^-(I,p). (A!) 

— OO — oo 

Since D ret (X, p) has only poles in the lower complex half-plane one can use contour inte- 
gration to evaluate the last integral. We assume that the self-energies can be neglected for 
large \p \, i.e. 

D Ict (X,p) -> — for \p \ -> oo. (A2) 
Po 

This allows us to write 

-oo \Ci+C 2 C 2 / \ / 

where C\ is a path along the real axis and C2 is an infinitely large half-circle surrounding 
counter clockwise the upper complex half-plane. 

For the effective spectral function we note that 1 — K can be written as (cf. (|4.2|)) 



with the ratio 



dr 

l-K = T— (A4) 
dpo 



p -£--ReZ^(X,p) 
r(X,p) = . (A5) 



Thus we get 



A= ¥ A ^-^ = kvrw^- (A6) 

Now the normalization of the effective spectral function can easily be shown: 

00 00 

^A= I—- - = 1 (A7) 

—00 —00 4 

where we have used that 

r — > ±00 for p — > ±00 . (A8) 

Again we have assumed that for large [_p 1 the self-energies can be neglected as compared to 
Po- 
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APPENDIX B: EQUATION OF MOTION FOR S 

In this Appendix we will derive (|4.8|) from ( ft.48 ). We have to show that 

[Po - ^ - ReE re \ S) - i [r, (po - £ - Re£ ret )S] 

- ~r\4 [p - ^ - ReS ret , S<] + i.4 [r, ( Po - f- - Re£ rct )S<] 

2 2m 2 2m 



vanishes. Using the definition 
[Po ~ 



V 2 - 1 

2m 2 J 



we find that (Bl) reduces to 

po - £ - Re£ ret 



r 



[r, s< 



(po-f-- ReS rct ) S< [log^- 
2m 

(po -|--I te g-) 5 < [ lo gW , - 



^-Re£ ret 1 
^ , 

1 1 , 



(Bl) 



(B2) 



with the ratio r defined in ( |A5| ). To finish our proof we only have to note that the last 
Poisson bracket in ( |B2j ) vanishes on account of the following general property of the Poisson 
bracket: 



[f 1 (g(X,p)),f 2 (g(X,p))] = 



(B3) 



for arbitrary functions g and arbitrary functions fi, f 2 which do not explicitly depend on X 
and p but only implicitly via g. 
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FIG. 1. (full) and (dashed) as functions of the energy po for the case of bosons. The 
width is chosen to be constant, T = 30 MeV. The temperature is T = 150 MeV and the kinetic 
energy = 100 MeV. See main text for more details. 



27 



